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1. S OLARS I M COMPUTER SIMULATION DEVELOPMENT 


SOLARSIM is a software package designed to predict the effect of 
certain electrical and mechanical imperfections of the SPS system on its 
performance. The following is a documentation of the SOLARSIM programs 
and SOLARSIM capabilities to quantify the spacetenna performance 
parameter values. The overall set up of SOLARSIM computer simulation is 
shown in Figure 1.1. As shown in the figure, the user input directs the 
central processor unit to nick the correct SOLARSIM subroutine from the 
subroutine package stored in the computer memory. Every operation done 
by the SOLARSIM is interactive in nature, i.e. , questions are asked and 
the user answers control the execution of the programs. As shown in the 
Figure 1.1, the SOLARSIM package can compute the following quantities: 

1. RMS Pointing error (PE) 

2. Tilt effects (TILT) 

3. MPTX code tracking loop performance (CDTL) 

4. MPTX carrier tracking loop performance (CRTL) 

5. Averaged power pattern (APP) 

6. Power transfer efficiency (PTE) 

There are two different types of inputs necessary for the operation 
of SOLARSIM routines, user inputs and the computer generated inputs. As 
mentioned earlier, the user inputs have the prompting, i.e., it is an 
interactive software, but the computer inputs are automatic in the sense 
that once the user's input specifies the SOLARSIM subroutine, the 
subroutine calls for the necessary computer generated input for- its 
execution. 

The SOLARSIM programs are configured such that all the user has to 
know is the set of inputs to the program he desires to execute. 
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Fig. Z.X VARIABLE AND FIXED SIZE SUBARRAY 
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Following is an attempt to make the user familiar with the assumptions 
and general requirements of the SOLARSIM package. 

2.0 SOLARSIM SPS SYSTEM DEFINITION AND COMPUTATIONAL CONSIDERATIONS 

2.1 Spacetenna 

Spacetenna is a stepped approximation of a circle of 11 cm in 
diameter having an ara of 0.76 km^. This area is subdivided into square 
subarrays totaling 101,552 N. The operation of the spacetenna being 

i 

retrodi recti ve, it needs a pilot beam originating at the center of the 
rectenna which is used to phase the downlink power beam to achieve 
retrodi recti Vi ty. The convention followed is as the direction 

V I 

of the pilot reference signal as seen from the spacetenna center. 

2.2 Current Taper 

For the control of the power lost in the si delobes of the antenna 
pattern, an excitation current taper of 10 dB (center to edge) is 
suggested. This taper is implemented by discretizing it into ten 
distinct levels and putting different numbers of amplifiers (klystron 
tubes) from level to level. It should be remembered that the ratings of 
all the amplifiers is the same all over the spacetenna. The ten levels 
when imposed over the circular spacetenna become ten circles with 
different diameters and having a different number of amplifiers per unit 
area from circle to circle. Figure 2.1 showi the spacet'^nna with ten 
power density circles, their radii and the number of klystrons in each 
circle. Even though the current rating of the klystron is the same, the 
current density is different in different power rings. The feed 
currents also have amplitude jitters associated with them. These two 
quantities are designated by: 

CURNT(I) Current densities in the power ring, 1=1,... ,10 
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SIGNAL(l) Amplitude jitter in the liil power ring, 1*1 10 

2.3 Subarrays 

The baseline system assumes only one power amplifier per 
subarray. Since there are different numbers of amplifiers present per 
unit area in different power density rings, the size of subarray is 
different from ring to ring. The subarrays are assumed to be square and 
have the power amplifier and the phase conjugation circuitry at its 
geometric center. The actual radiating elements are the slots in the 
waveguides. Figure 2.1 shows the size and number of subarrays in each 
power ring. One alternate way is to have a constant sized square 
subarray throughout the spacetenna. Such an arrangement will have a 
different number of power amplifier per subarray from r-ing to ring. 

This arrangement is also shown in Figure 2.1. The following convention 
is adapted: 

N{I) The number of radiating slots in i— power ring 


1 = 1,... ,10 


M(I) The number of radiating slots per subarray in i^ 

power ring, 1=1,... ,10) 

2*4 Tilts on Subarrays 

There are basically one type of tilt and two types of jitters 
associated with the subarrays which are, mechanical tilting of the 
subarrays with the associated jitters and the radiating and transmiting 
element location jitters. The mechanical tilittng of the subarrays can 
be considered to have two components, i.e, the x-component and the y- 
component the mean values of which gives the mechanical tilt of the 
spacetenna with respect to the perfectly pointing spacetenna. The 
radiating and transmitting element location jitters actually have three 
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Fig. 2.2. EFFECTS OF SUBARRAY MEAN TILTS AND JITTER 
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components: two components in the plane perpendicular to the line 
forming that particular eleuent and the rectenna center and one 
component along that line. We will neglect the jitter i )mponents in the 
plane mentioned above and consider only the jitter component along the 
line joining the element and the rectenna center as the location jitter 
of that particular element. The convention followed is given in the 
following table: 


Mean Mechanical Tilts of the 
Subarrays 

Jitters on Mechanical Tilts 

Location Jitter on the Pilot 
Receiving Element 

Location Jitter on the Radiating 
El ement 


SIGMA(I), I»l,2,...,10 
SIGPSI(I), 1*1,2 10 

SIGPHI(I), '1*1,2,. ..,10 


It should be noted that as well as minutes while 

the location jitters will be specified in terms of It of the wvelength of 
the power wave. Figures 2.2 and 2.3 show these tilts and jitters. 

2.5 Phase Control System 

Since the operation of the spacetenna is retrodirective, it needs a 
constant phase reference throughout the antenna for conjugator to 
function properly. This constant phase is supplied to the conjugators 
by the use of MSRTS in the form of a tree. The master oscillator 
situated at the physical center of the spacetenna locks onto the phase 
of the incoming pilot signal. This phase is transmitted to 16 first 
level slave oscillators, each of these sends the phase to 16 second 
level slave oscillators. These third level slaves in turn sene the 
phase to 25 fourth level slaves. Figure 2.4 shows the phase 
distribution tree. 
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As the master phase passes from level to level of the phase 
distribution tree structure it gets corrupted by the phase noise added 
by the oscillators and the related equipment at each level. The 
variance of the accumulated noise is designated as SI6MAB(I), at the end 
of the distribution tree. 

2.6 Inputs for the Double Integrator Subroutine Supplied by the 
Computer Library "Math PacP~ 

The description of these inputs is necessary for the power transfer 
efficiency program only because a double integration of the averaged 
power pattern is necessary to obtain the total power radiated by the 
spacetenna and the power received by the rectenna. The following 
formula will illustrate the Inputs necessary for the program: 


Power Received 
on Area 'A' 


e a 
u ’^u 

/ / APP(e,^) sin eded^ 


where APP(0,4») is the averaged power pattern produced by the 
spacetenna. Where 

®L*®u* lower and upper limits of the variable 9 

describing the area A. 

The lower and upper limits of the variable <f 
describing the area A. 

Note: To obtain the power received by the rectenna 9 l=0, 

9y = 0.477', = 0 and = 360°. 

The double integral subroutine approximates the above double 
integral by the following double sum 


9 

u 


/ 


i 

P 

/ APP(9,^)sin 6d9di> = 

♦l 


I 


^ % 


k=l t=l i=l j=l 


l 1 8k' 
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where 

Are the number of intervals on the 6 axis and ^ axis 
respectively, on the specified e and ♦ ranges. 

Ng,N^: Number of points in each of Ig and respectively in 

the e and ranges. 


^6i‘ \i “^tj generated by the subroutine and does not have 
to be supplied. 

Note: A reasonable input for Iq» 1^ would be 3,3 and for Ng,N^ wo.uld be 
4,4. 

3. USAGE OF SOLARSIM SUBROUTINES 

SOLARSIM subroutine package is created such that user encounters 
minimum of trouble to run it. After a few prel i mi nary 'commands from the 
user, the entire operation becomes automatic stopping and asking 
questions where user input is necessary. User is prompted for all the 
inputs necessary for that particular program; thus, virtually 
eliminating all user generated errors in data feeding. The preliminary 
commands preparing the SOLARSIM subroutine package for user selection of 
program are necessary only once, making selecting and running of 
subsequent programs easy. The preliminary commands necessary are shown 
in the following computer printout. As seen on that sheet the last 
command is 0XQT LINCOM. SELECT. This executes the 'select* subroutine 
which shows that there are six choices possible. They are: Averaged 

Power Pattern, Power Transfer Efficiency, Pointing Error, Tilt Affects, 
Carrier Tracking Loop and the Code Tracking Loop. @ ADD. LINCOM. 
Abbreviated File Name is the command necessary to access the required 
subroutine from the magnetic tape which was loaded on the tape drive 
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3.1 SOLARSIM Subroutine Pointing Error (PE) 

Purpose : SOLARSIM subroutine POINTING ERROR evaluates the effect of 
phase error introduced by the phase distribution tree (at va^'ious 
levels of the tree) onto the pointing error. It also allows 
variations in number of levels and number of branches per node at 
each level of the phase distribution tree. This subroutine may 
also be used for determining effects of other phase error sources 
on the pointing error. An example of this would be the phase error 
introduced by the SPS transponder circuitry can be counted as the 
pahse error due to the last level of the phase distribution tree. 
Access Command: 13 ADD LINCOM.PE 

Inputs : The only necessary inputs are: 

Number of levels in the phase distribution tree (PDT) 

Number of power density levels 
Number of branches per node at each level of PDT 
Out put : The following page shows a sample run of the program. 
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3.2 SQLARSIM Subroutine TILT 

Purpose ; SOLARSIM subroutine TILT evaluate the effect of various tilts 
of the subarrays, the location jitters 3f the power radiating and 
pilot receiving elements and finally the phase distribution system 
errors on the gain of the spacetenna. 


Access Comnand: 

0 ADD LINCOM.TILT 

Inputs: 


CURNT(I) 

I ^ 1,...,10 

N(I) 

I = 1 10 

M(I) 

I = 1 10 

SIGPSI(I) 

I » 1,...,10 

SIGPHI(I) 

I = 1,...,10 

YMEAN(I) 

I = 1 , . . . jlO 

XMEAN MAX 

Maximum XMEAN tilt in minutes 

SIGMA(I) 

I lif2,.^.,10 


Note that this program does need the integrator subroutine to generate 
the weights. One needs to know the above inputs only if values 
different from baseline values are to be fed to the program. 

Outputs: The following pages show a sample run of the program. 
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3.3 SOLARSIM SubroMtine Power Pattern (AP?) 

Purpose ; The purpose of this subroutine is to study the effects of 
parameters like the subarray tilt>j, the phase jitters, ';urrent 
amplitude jitter, etc. on the power pattern produced by the 
spacetenna. 

Access Command; @ ADD LINCOM.APP 

Inputs ; 


N(I) 

I = 

1,...,10 

M(I) 

I » 

1,...,10 

NA(I) 

I = 

1,...,10 

R(I) 

I = 

1,...,10 

SIGMAB(I) 

I * 

1,...,10 

SIGPSI(I) 

I = 

1,...,10 

SIGPHI(I) 

I = 

1 10 

SIGMAI(I) 

I = 

1 10 

CURNT(I) 

I = 

1,... ,10 

XMEAN(I) 

I = 

1,...,10 

SIGMA(I) 

I = 

1 10 

PHI 

The 

value of 


generated. 

NTERM 

The 

number of 


be evaluated. 

STEP 

The step size 


Note that the parameter values N(I), M(I), and R(I) are required to be 
fed into the program only if these values are different from the 
baseline values. All the necessary baseline values are already present 
in the program. Of course, NTERM and STEP will be required to be fed in 
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by the user. 

Output ; The following pages show a sample run of the subroutine APP. 
In this particular run the total rms phase error is assumed to be 
10° (input) while all the other inputs were held at zero. It 
should be noted that the program also produces the power 
output. This output is true for the SPS at 36,000 km from the 
surface of the Earth. 


- 21 - 


P. ^ 

caLinL^m 





In 

Oi rH m 

r> fi o 

mp o 
f? O 
01 o 

3 

V 

r* f> f-< n 

IM m n o 

•r f> Vfl ' 
fS P o < 
r«l n 
m ^ 

40 « 


O p GA 
O O P I 
P o o < 
p e o f 
? • r f 


o o o 
o o o 
o o o 
o o o 


V 

•<-< <o 
rM o CO o 
« cv ^ ^ 

r-4 «M m 


o o o 

OOP 

OOP 

poo 


Ip 

00 ri f** fo 
«-i cn (o o 

P lO iO 4- 
#H r4 o 


o o o 
o o o 
OOP 
poo 


P ot 


40 p o ro 
rj rt to o 
#n o o 4> 

f-4 «r r4 CD 
f-< 40 

to 

2 

r4 • 

f4 rr ^ n 
O 40 o 
O P fX 4- 
jH to Hi ro 
to o 

ifi (N 

f4 m 

V 

o ^ p m 
IN O 40 o 
r* N «s 4- 
lO H H I 


O O O 

o o o 
o o o 

O CD O 


» fT1 O O p 
Ml) P P 7 
» O P O to 

> «r o o p 

j» • • fX 

1 I 

> 4- H 

> "T O P P 
I H O p 4^ 

r> o o H 

I p p a to 

* • • p 

: I 

POPP 
O P P f 
o o o rf 

1 CD o rr* 

"i 

O 

4' 'H 
m o o p 
fr o p 4 
fr O O 0 
40 p p 4t 

> • • o 


O P P 
OOP 
POO 
OOP 


4 

POOP 

p o o 4 
p o o «r 
4D o o p 

> • *P 


POO 

OOP 

OOP 

pop 


o o o 

40 O o 

•r o o fM 
rr o p m 


^ • • #f>opm 

M I * • • (TJ 

to <X IX 

o o 

O i P I P 

OH040I0 0 004 H 

ooa^«oOQO<0Oop 
Op«rP4'0004^P04 
O pHtoOOOHOOto 

O f> Hf**toOOIX 

o r> r) ^ 


to IX 

1 • 


p o p ♦ 
o o o to o o 0 

o o a ry o o 4 
p o o H p o 4 
» * *0000 


. { *=* 

: p H p in 
MX to to o 
; 40 p to 4 
f to H n 
p 

to H 


0 4o4*^*^0004 

Otop.toOOOPO< 

o(nrxto4‘pooo4 

O to P P O O P I 

O H P* k • • i-f < 

gas: 1 


> P O O P O O ( 
• • rf o O r 
• ' • •! 
i r 

I ■ ■ 

^ ^ ^ I i N 
M M H ^ 


toHH'^-'Mpouact 

4ta:s;zrtftocotoU» 




1 

j 



ENTER THE STEP SIZE IN MINUTES: 

> 1.0 



Purpose ; The purpose of SOLARSIM subroutine power transfer efficiency 
is to study the effects of perturbations, aiechanical or otherwise 
on the power transfer efficiency of the spacetenna. 

Access Command: 0 ADD LINCOM.PTE 

Inputs ; If the program ,1s run for the baseline quantities then the 
following quantities are not necessary but if they are to be 
different from baseline numbers one needs it for the input of the 
program. 

N(I) I = 1 N 

M(I) I = 

NA(I) I = 

R(I) I = 1 N 

SIGMAB(I) I = 1 N 

SIGPSI(I) I = 1 N 

SIGPHI(I) I = 1 N 

SIGMAI(I) I = 1 N 

CURNT(I) 1 = 

XMEAN(I) I = 

SIGMA(I) I = 1 N 

The pilot incidence angles as seen from the 
spacetenna. 

Xa.Xb Lower and upper limits of theta integration in 

mi nutes. 

Lower and upper limits of phi integration in, 
degrees. 


Ya.Yb 






Nx,Ny 


Number of points to be used per interval in the 
Gauss Quadrature approximatin of the integral 
Ix.Iy Number of intervals in the theta and phi range to be 

considered for the integration. 

Output : The following pages indicate the computer run for the Power 
Transfer Efficiency Program. In this particular run the location 
jitters on the radiating and receiving elements as well as the 
current amplitude jitters are constants for the spacetenna, i.e. , 
they do not change from power ring to power ring. In such a case 
the program has the capability of looping back and changing the 
initial conditions of the parameters. As seen from the printout, 
the program is run for different initial conditions for the 
parameters. Any one of the parameters is allowed to vary producing 
a set of values for the efficiency holding all other parameters to 
constant values. 
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3.5 SOLARSIM Subroutine Carrier Tracking Loop (CRTL) 

Purpose : The purpose of SOIARSIM subroutine carrier tracking loop is to 
study the effect of the random various noises and various system 
parameters such as the filter bandwidth, loop bandwidth and the 
chip rate on the carrier tracking loop phase error. 

Access Command: 0 ADD LINCOM.CRTL 

Inputs : The question and answer session generates the following Inputs: 
f3 “ RF filter 3 dB bandwidth 

6 = Half 3 dB notch frequency normalized to the chip rate 

Y = The front end filter null in per unit value 

3 = PN apparent tracking offset normalized to chip time 

M = Code repetition rate 

Rj, = Chip rate 

kj = Noncoherent interference coupling coefficient 

k2 = Coherent interference coupling coefficient 

It should be noted that the f3 ^(0 value should be 0.0 when there is no 
RF filtering before the front end filter. This is provided in the 
program to account for the filtering effect induced by the waveguides in 
the subarray. Also the parameter 3 is not used in the program because 
it is assumed that there is 10% loss of pwoer due to the code tracking 
error. The results are rather pessimistic because the code tracking 
loop can track the code very well. 

Output : The computer run of SOLARSIM subroutine CRTL is shown in the 
following pages. 
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Purpose; The main purpose of the SOLARS IM subroutine code tracking loop 
is to study the behavior of the code tracking loop under various 
noises. The subroutine may be used to produce design values for 
the code tracking loop such as the arm filter bandwidth or the code 
loop bandwidth* etc. 

Access Command ; 0 ADD LINCOM.CDTL 
Inputs ; This program has only three inputs; 

Bl * The loop bandwidth in Hz 

^3 <fi “ One-half 3 dB bandwidth in kHz of the bandpass filter 
in the code loop arm. 

“ The dither frequency in kHz. 

Ncte; F 3 (jB is half of the 3 cB bandwidth and not the total bandwidth. 
If this instruction is not followed, erroneous results will be 
generated. 

Output ; The following pages show a sample run of the SOLARSIM 
subroutine CDTL. 










